Competition of resonant and nonresonant paths in resonance-enhanced two-photon 
single ionization of He by an ultrashort extreme-ultraviolet pulse 



(N 

o 

(N 



O 
Ph. 



> 

On 

m 



3 



Kcnichi L. Ishikawa x '0 and Kiyoshi Ueda 2 

1 Photon Science Center, Graduate School of Engineering, 
The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan 
2 Institute of Multidisciplinary Research for Advanced Materials, Tohoku University, Sendai 980-8577, Japan 

(Dated: January 19, 2012) 

We theoretically study the pulse-width dependence of the photoelectron angular distribution 
(PAD) from the resonance-enhanced two-photon single ionization of He by femtosecond (< 20 fs) 
extreme-ultraviolet pulses, based on the time-dependent perturbation theory and simulations with 
the full time-dependent Schrodinger equation. In particular, we focus on the competition between 
resonant and nonresonant ionization paths, which leads to the relative phase 8 between the S and D 
wave packets distinct from the corresponding scattering phase shift difference. When the spectrally 
broadened pulse is resonant with an excited level, the competition varies with pulse width, and, 
therefore, 5 and the PAD also change with it. On the other hand, when the Rydberg manifold is 
excited, 8 and the PAD do not much vary with the pulse width, except for the very short-pulse 
regime. 

PACS numbers: 32.80.Rm, 32.80.Fb, 41.60.Cr, 42.65. Ky 



Multiphoton ionization of atoms has consistently been 
receiving a great deal of attention for decades (see e.g.[l|- 
ll|). The advent of intense extreme- ultraviolet (EUV) 
sources such as high-harmonic generation (HHG) and 
free-electron lasers (FEL) has enabled two-photon ion- 
ization (TP I) of sp ecies with a deep ionization potential 
such as He (l2l - flq and N2 [l7|- Upon photoionization, 
the continuum electron wave packet is emitted, which 
is a superposition of different partial waves, each with 
its own orbital angular momentum, intensity, and phase. 
Photoelectron angular distribution (PAD), nowadays ex- 
tensively studied by the velocity map imaging technique 
(see e.g. HIE!), contains information on the interfer- 
ence of these different partial waves. 

In this Letter, we theoretically study the pulse-width- 
dependence of the PAD from two-photon single ioniza- 
tion of He by femtosecond (fs) EUV pulses. Especially, 
we focus on situations where the pulse is closely reso- 
nant with an excited level, i.e., resonance-enhanced TPI. 
We have chosen He as a target atom for the following 
reasons: first, its single-electron excitation energies, e.g ., 
21.218 eV for ls2p 1 P and 23.087 cV for ls3p 1 P |2fl|, 
coincide with the 13th and 15th harmonic photon ener- 
gies of a TkSapphire laser, respectively, and also with 
the typical wavelength range of EUV FELs such as the 
Spring-8 Compact SASE Source (SCSSjEH, the Free- 
electron LASer at Hamburg (FLASH) J^jfTand FERMI 
[23| . Second, its simple electronic structure allows for ex- 
act time-dependent numerical analysis [24 - 28 1 , in great 
contrast to alkali atoms. 

In the case of resonance-enhanced TPI, the resonant 
ionization path via resonant levels and the nonresonant 
path via nonresonant intermediate levels coexist 0] ■ Our 
results show that in the few fs regime, the competition 
between the two paths can be controlled by changing 



the pulse width when the pulse is resonant with a single 
excited level. The relative phase 6 between the differ- 
ent partial waves (S and D for He) would be just the 
scattering phase shift difference for nonresonant TPI and 
resonant (1 + 1') TPI [lj||. For the present case, on the 
other hand, the ionization-path competition give rise to 
an additional contribution S ex . This contribution clearly 
manifests itself in the PAD, and both the PAD and <5 
vary with the pulse width. When the pulse becomes 
so short that its spectrum gets broader than the level 
spacing and resonant with multiple levels, especially the 
Rydberg manifold, 5 is still different from the scatter- 
ing phase shift difference, but does not vary with T. We 
further explore how the chaotic nature of FEL radiation 
2^33 1 affects the PAD. Our analysis using the partial- 
coherence method [HI indicates that the PAD is between 
those corresponding to the coherence time and the mean 
pulse duration. 

Let us first do a simple analysis on how the relative 
importance of the resonant and nonresonant paths de- 
pends on pulse width, based on the second-order time- 
dependent perturbation theory within the common ro- 
tating wave approximation. The dynamic Stark effect is 
negligible for pulse parameters used in the present study. 
We consider the process where a laser pulse with a central 
frequency ui and a pulse envelope f(t), linearly polarized 
in the z direction, promotes an atomic electron from an 
initial state \i) to a final continuum state |/) through 
two-photon absorption. The complex amplitude Cf of 
the final state after the pulse in the interaction picture 
can be written as, 
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where fi m n denotes the dipole transition matrix ele- 



2 



merit between states m and n, A m = uj m — (uii + uS), 
Af = ujf — (uj rn + lu) with ui m being the energy eigen- 
value of state m, and the sum runs over all the inter- 
mediate bound and continuum states m. Although a 
rectangular pulse is often assumed in previous work 
we take, as a more realistic choice, a Gaussian profile 
/(*) = E e~ t2 / 2T \ with E a and T being the field ampli- 
tude and the pulse width, respectively. For the case of 
ujf = cji + 2uj, in particular, one can perform the integrals 
in Eq. (JlJ analytically to obtain a physically transparent 
expression: 



c f = ttE^T 2 t*fr. 
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where F(x) denotes Dawson's integral [35[, which tends 
to x near the origin and l/2x for x — > oo. Only resonant 
states within the spectral width of the pulse contribute to 
the first term, corresponding to the resonant path. On 
the other hand, the asymptotic behavior of F(x) sug- 
gests that all the intermediate states except for the exact 
resonance (A m = 0) participate in the second term, as 
expected for nonresonant paths. While either term domi- 
nates for a relatively long pulse (ps and ns) , we can expect 
that the two terms are comparative for ultrashort (~ a 
few fs) pulses and that their relative importance, which 
may be expressed as argc/, varies with T. In such a sit- 
uation, the amplitude ratio cs/cd between the final S 
and D continuum states is complex, since the branching 
ratio fism/ fJ-Dm of the transitions from the intermediate 
P states m to each state depends on m. While the actual 
outgoing wave packets involve the contribution from the 
final states with w/ ^ u>i + 2uj, it is instructive to write 
argcs/co using Eq. ^ as follows: 
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With a f (f = S,D) = (fif r fl ri ) 1 J2m.(^r) HfmHmi/A m . 

Here we have assumed that only one intermediate state 
r is resonant with the pulse and that F(A m T) w 
(2A m T) _1 for all the other intermediate states. Hence, 
the competition between the resonant and nonresonant 
paths affects the interference between the outgoing S and 
D wave packets and manifests itself in the photoelectron 
angular distribution. 

The photoelectron angular distribution is given by [36| , 
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where a is the total cross section, 6 is the angle between 
the laser polarization and the electron velocity vector, 
and /3 2 and ,84 are the anisotropy parameters associated 
with the second- and fourth-order Legcndre polynomials, 
respectively. The interference of the S and D wave pack- 
ets produces a photoelectron angular distribution propor- 
tional to \\cs\e lSo Y 0() - \c D \e ,S2 Y 20 \ 2 , with Si being the 



phase of the partial wave, or the apparent phase shift. 
Then, the anisotropy parameters can be described by, 
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where W — \cs/cd\ and 5 = So — S2 j37|. The apparent 
phase shift difference S = S sc +5 ex consists of a part 5 SC in- 
trinsic to the continuum eigen wave-functions (scattering 
phase shift difference) . which has previously been studied 
both theoretically (38l - l40l | and experimentally [3], and 
the extra contribution S ex = arg cs/cd from the competi- 
tion of the two paths. This situation present a contrast to 
the case of the photoionization from photo-excited states 
[HI], where the nonresonant path is absent and only S sc 
is present (S = S 8C ). 

We now verify this qualitative idea, using direct numer- 
ical solution of the full-dimensional two-electron time- 
dependent Schrodinger equation (TDSE) (28^ : 



^ $ (y,*) =[go+gf(t)]$(riir2>t)i 

with the atomic and interaction Hamiltonian, 
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H T (t) = (z 1 + z 2 )f(t)smut. 
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We solve Eq . ([6]) using the time-dependent close-coupling 
method 24|- 2^ . The numerically obtained excitation en- 
ergies for the \s2p 1 P and ls3p 1 P states are 21.220 and 
23.086 eV, respectively, in fair agreement with the exper- 
imental values (21.218 and 23.087 eV [2(J, respectively). 
Sufficiently after the pulse has ended, we calculate /3 2 
and /?4 by integrating the ionized part of |<f>(ri, r 2 )| 2 over 
Ti , r 2 , # 2 , 0i , 4>2 , from which one obtains W and S by solv- 
ing Eq. (j5|) . We use the values of S sc from to calculate 
Sex = S — S sc . The calculation has been done for a Gaus- 
sian pulse envelope with a peak intensity of 10 11 W/cm 2 , 
at which we have confirmed that the interaction is still 
in the perturbative regime. 

The pulse-width dependence of 6 and W for hw = 
21.2 eV and hu = 21.3 eV close to the \s2p reso- 
nance (21.218 cV) is shown in Fig. [T] The calculations 
have been done at different values of full-width-at-half- 
maximum (FWHM) pulse width T 1/2 = 2Vm~2T be- 
tween 500 as and 21 fs. As expected, both 5 and W 
substantially changes with pulse width, especially when 
the pulse is shorter than 10 fs. Accordingly, the PAD also 
varies as shown in Fig. [2j One finds that the distribution 
to the direction perpendicular to the laser polarization, 
i.e., 8 sa 90°, 270° decreases as the pulse is shortened. 
This can be understood as follows: roughly speaking, 5 
changes from ~ ? to ~ n as T]/ 2 varies from 21 fs to 
500 as. Thus, c s /cd is approximately real and negative 
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FIG. 1: (Color online) The pulse- width dependence of the 
TDSE-derived apparent phase shift difference (relative phase) 
S (left axis) and W = |c s /cd| (right axis) for Titjj — 21.2 eV and 
21.3 eV. The thin horizontal line denotes the value (1.511 
[4(J) of the intrinsic scattering phase-shift difference 5 SC for 
hu = 21.2 eV. That for hui = 21.3 eV is 1.491 Solid lines 
are the results of fitting using Eq. (3). 
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FIG. 2: (Color online) The pulse- width dependence of the 
photoelectron angular distribution for hui = 21.2 eV. 



in the short-pulse limit, which leads to the cancellation 
between Yqq(9, <p) and ¥20(8, cp) around 8 = \. As stated 
earlier, strictly speaking, Eq. §3§ is applicable only to 
LUf = uji + 2ui, and the actual PAD involves integration 
over ujf. Nevertheless, the results in Fig. Q] can well be 
described by Eq. ([3]) (solid lines in Fig. [T]), except for 
S in the ultrashort pulse regime Ty 2 < 1 fs, where the 
spectrum becomes broader than the level spacing. 

With increasing pulse duration, 5 approaches the scat- 
tering phase shift difference S sc , and the PAD changes 
only slowly with Tyz (Fig. [2]). When the pulse is reso- 
nant (A r T <C 1) and sufficiently long (T 3> as,ao) at 



FIG. 3: (Color online) The spectral-width dependence of the 
TDSE-derived apparent phase shift difference (relative phase) 
S for hui — 21.2 eV. The thin horizontal line denotes the value 
(1.511 [i(|) of 8 ac . Thin dashed line plots the asymptotic 
behavior Eq. ([9j). 



the same time, assuming that the resonant excitation is 
not saturated, one can approximate the extra phase shift 
as 



(a D - a s )/\/TrT, 



(9) 



hence, it is proportional to the spectral width, which can 
be confirmed in Fig. [3] 

On the other hand, if we plot S as a function of spec- 
tral width (Fig. [3]) , S tends to an asymptotic value in the 
wide-spectrum, i.e., short-pulse limit. Correspondingly, 
the PAD does not change much with the pulse width 
for Ti/2 < lfs (Fig. [2]). This is because the pulse be- 
comes resonant with multiple levels; the spacing between 
the ls2p and ls3p is 1.9 eV. Indeed, one can show that 
when many neighboring states are resonantly excited by 
the pulse, the extra phase shift difference S ex does not 
much depend on the pulse duration. This especially ap- 
plies when the photon energy lies in the Rydberg man- 
ifold, and exceeds the ionization potential (24.59 eV), 
i.e., in the case of above-threshold two-photon ioniza- 
tion. In Fig. |4] we compare the pulse- width dependence 
of Sex for different values of hui. While 6 ps S sc [8 ex ~ 0) 
for nonresonant pulses (Hui = 20.3 cV, Ti/ 2 > 3.5 fs), 
when the pulse is close to resonance with an excited 
level (Hui = 21.2,21.3, and 23.0eV), S ex rapidly changes 
with Ty%. On the contrary, S ex is nearly constant for 
Hu = 24.3,24.6, and 25.0 eV. At T 1/2 < 2 fs, the spec- 
trum is so broad that S ex restarts to change slightly. One 
also sees that the transition across the ionization poten- 
tial is smooth. It should be pointed out that the extra 
phase shift difference due to free-free transitions plays 
a significant role in the recently observed time delay in 
photocmission by attosecond EUV pulses 
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Coherent EUV pulses considered so far may be gen- 
erated by HHG sources and HHG-seeded FELs 0. As 
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FIG. 4: (Color online) Pulse-width dependence of S ex for dif- 
ferent values of photon energy. The resonant excitation en- 
ergies for ls2p 1 P, ls3p 1 P, ls4p_P, ls5p x P levels are 21.218, 
23.087, 23.742, and 24.046 eV M. 



is well known, on the other hand, the temporal pulse 
shapes of FEL operating in the self-amplified sponta- 
neous emission (SASE) mode fluctuate from shot to shot 
In order to investigate the effects of the chaotic 
nature, we have performed numerical experiments for 
EUV pulses randomly generated by the partial-coherence 
method j34|. The obtained values of W and S for sev- 
eral combinations of coherence time (CT) and mean pulse 
width (MPW), and the corresponding PAD are shown in 
TableUand Fig. 03 respectively. We can see from the table 
that 6 and W take values between those corresponding 
to the CT and the MPW in most cases. Also, the PAD 
is of a shape between those for the CT and the MPW. 
The detailed mechanism underlying these somewhat em- 
pirical findings will require further investigation taking 
photon statistics into account. 

In summary, due to the competition between the reso- 
nant and nonrcsonant ionization paths, the relative phase 
5 between the S and D photoelectron wave packets from 
the resonance-enhanced two-photon ionization of He by 
fs EUV pulses is different from the scattering phase 
shift difference S sc which would be expected for single- 
photon ionization 18[ and nonrcsonant two-photon ion- 
ization, and rapidly changes with the pulse width when 
the pulse is resonant with an intermediate excited state 
and 2fs < T\/ 2 < 10 fs. Accordingly, the photoelectron 
angular distribution varies with 2\/2 as well. Also, 5 ex 
is finite but constant independent of T1/2 when the Ryd- 
berg manifold is excited . Hence, the control of the com- 
petition between the resonant and nonresonant paths in 
He by pulse width is a unique feature of a-few-fs EUV 
pulses. The PAD is affected by the chaotic nature of 
SASE FEL pulses, and takes a shape between those corre- 
sponding to the coherence time and the mean pulse width 
of the pulses. The results of the present study stress 
the importance of the account of the nonresonant paths 



TABLE I: W and S calculated for chaotic pulses generated 
by the partial-coherence method for several pairs of coherence 
time (CT) and mean pulse width (MPW). The average values 
and standard deviation errors of W and 8 by 48 runs are listed. 
The rows with the same CT and MPW values are for fully 
coherent pulses as in Figs. [T]and[4] 



hw (eV) CT (fs) MPW (fs) 
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21.2 


2 
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1.31 ±0.08 


2.03 ± 0.04 
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1 42 + OQ 


1 91 ± 04 
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1.49 ±0.09 


1.87 ±0.03 




2 
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1.92 
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0.645 ± 0.061 


2.14 ±0.08 
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0.890 ± 0.058 


2.15 ±0.06 
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2.45 




7 


7 


1.08 


2.10 



180 




270 



FIG. 5: (Color online) Photoelectron angular distribution by 
chaotic pulses for hw = 21.2 eV and three pairs of (CT)- 
(MPW) indicated in the legend. 



in the interpretation of resonant two-photon, single- or 
two-color, ionization experiments by state-of-the-art ul- 
trashort EUV sources. 
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